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> : Abstract 

^ ! We discuss the possible impact of astrophysical neutrino flux measurements at neutrino 

O [ telescopes on the neutrino oscillation program of reactor experiments and neutrino beams. 

^ I We consider neutrino fluxes from neutron sources, muon damped sources, and pion sources, 

■~>^ ■ where we parameterize the input from these sources in terms of the flux ratio R = 0^/(0e + 

p ■ (pr) which can be extracted from the muon track to shower ratio in a neutrino telescope. 

JL , While it is difficult to obtain any information from this ratio alone, we demonstrate that 

D ! the dependence on the oscillation parameters is very complementary to the one of reactor 

'^ ' experiments and neutrino beams. We find that for large values of sin^ 26'i3, a measurement 

of R with a precision of about 20% or better may not only improve the measurement of the 

^ ■ leptonic CP phase, but also help the determination of the mass hierarchy. In some cases, 

early information on 6cp may even be obtained from Double Chooz and an astrophysical flux 

alone without the help of superbeams. For small values of sin^ 26'i3, we flnd that using the 

information from an astrophysical neutrino flux could eliminate the octant degeneracy better 

than reactor experiments and beams alone. Finally, we demonstrate that implementing an 

additional observable based on the electromagnetic to hadronic shower ratio at a neutrino 

telescope (such as at higher energies) could be especially beneflcial for pion beam sources. 
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1 Introduction 

In the coming ten years, many experiments with a terrestrial neutrino source will provide 
information on the neutrino oscillation parameters. In particular, the parameters sin^2^i3, 
(5cp, and the neutrino mass hierarchy are very interesting, since we know very little about 
them except from an upper bound on sin^2^i3 [1]. However, even if sin^ 2^13 is large, mea- 
suring leptonic CP violation and the neutrino mass hierarchy with long-baseline and reactor 
experiments will be extremely challenging on that time scale because of the low statistics 
(see, e.g., Ref. [2]). Therefore, any complementary information could help to improve the 
measurement. In this study, we discuss the potential complementarity of three sources 
of information: Neutrino beams, reactor experiments, and a possible flavor composition 
measurement of an astrophysical neutrino flux. 

Considered sources 

Neutrino beams, such as MINOS [3], T2K [4], or NOz/A [5], are sensitive to sin^ 26*13, ^cp; 
and the neutrino mass hierarchy via electron neutrino appearance. However, this ability to 
access all of the parameters leads to correlations and degeneracies [6-9] when it comes to the 
extraction of the individual parameters. In addition, the next generation of "narrow-band 
beams" operated at the oscillation maximum will be mainly sensitive to the CP-odd part 
of the appearance probability. Reactor experiments [10], such as Double Chooz [11], will 
be very complementary to provide "clean" information on sin^ 2^13 only [12,13]. However, 
these experiments cannot access (5cp and the mass hierarchy at all, and the statistics of the 
first generation of experiments will be moderate. 

As another potential source of information, a detection of astrophysical neutrinos at neutrino 
telescopes [14-17] with a well-predicted flavor composition at the source (such as from 
neutron or pion decays) could provide additional knowledge on the mixing parameters [18- 
20]. The existence of such astrophysical neutrinos is not yet proven, but the detection of very 
high energy cosmic rays suggests that cosmic accelerators exist which should also produce 
high energy neutrinos, such as by pion decays. There are many potential candidates for 
neutrino sources, such as gamma ray bursts, active galactic nuclei, or starburst galaxies, 
where the latter source is not affected by the Waxmann-Bahcall bound [21]. Since we are 
not interesting in the type of source, but the initial flavor composition in this study, we 
will follow a more abstract approach and classify potential sources by their initial flavor 
composition below. If the distance to the source is long enough or the production region 
is large enough, neutrino oscillations on the way from the source to the Earth will average 
out. The predicted flavor composition at the Earth then depends on the mixing parameters 
including (5cp and the CP-even part of the mixing only (see, e.g., Ref. [22]), whereas it is 
not sensitive to the mass hierarchy. However, compared to the terrestrial neutrino sources, 
there will be no spectral information, which means that, given the low statistics, it will be 
difficult to extract some information on the individual mixing angles. 

For the astrophysical sources, we consider neutrino production from neutron decays with 
the flavor ratio (1:0:0) at the source [23,24] ("Neutron beam source"), from pion decays 
with the flavor ratio (1:2:0) at the source ("Pion beam source"), and from muon damped 



pion decays with a flavor ratio (0:1:0) at the source [25,26] ("Muon damped source"). In 
the last case, the muons are absorbed before they can decay, which may also occur for pion 
beams at high energies [26]. The neutrino fluxes expected from these sources typically scale 
as dN/dE oc -E^^, where the typical upper neutrino energies can be very different.^ For 
example, for gamma ray bursts, the typical neutrino energies depend very much on the 
production region (and process) and can be between GeV and EeV. For a neutron beam 
produced by photo-dissociation of heavy nuclei, the typical energies are more limited to the 
TeV range. The sources can be either steady (such as for a neutron beam produced by photo- 
dissociation of heavy nuclei, a pion beam from active galactic nuclei, or a diffuse background 
flux), or time-dependent on a short timescale (such as for a flux from a gamma ray burst). 
Depending on neutrino energies and steadiness of the source, the background especially from 
atmospheric neutrinos can be problematic for the event extraction especially below about 
1 PeV. For more discussion on background discrimination approaches, see Appendix 1X1 

Observables and statistics for astrophysical neutrinos 

As far as the detection of the astrophysical neutrinos is concerned, we do not distinguish 
between neutrinos and antineutrinos, i.e., the neutrino and antineutrino fluxes are simply 
added. A neutrino telescope can identify muons by their tracks (Cherenkov light with con- 
tinuous loss of energy). However, electron and tau neutrino events are harder to disentangle 
because they both produce showers of particles with a larger threshold. Depending on the 
neutrino energy, the electromagnetic showers (from electrons) and hadronic showers (from 
taus) may be disentangled by their muon content, as well as the tau track may be measured, 
which means that there could be a possibility to disentangle electron and tau neutrino neu- 
trino events (see, e.g., Refs. [27,28]). An additional process to identify Ve^ is the Glashow 
resonant process z/g + e" — > W~ at 6.3 PeV, which can also be used for oscillation parameter 
measurements [20]. For most of this study, we follow Refs. [18,19] (see also Ref. [29]) and 
use R = (j)fj,/{(j)e + (j)r) as observable, which means that we make the conservative assumption 
that electron and tau events cannot be disentangled. This observable can be extracted from 
the ratio of muon tracks to showers [28] , but note that there is an additional hadronic shower 
background from neutral current events for all flavors which needs to be subtracted. We 
assume that we have a precision measurement of R with an effective relative error exploiting 
all available information and containing all systematics and backgrounds. For example, the 
information from different energies and event types may be used to reduce the systemati- 
cal errors on that quantity (see, e.g., Refs. [26,30]). Especially using different energies is 
plausible for neutrino telescopes, since we assume that neutrino oscillations are averaged 
out (compared to the beams and reactor experiments, there will be no energy dependence 
of the oscillations). Finally, we assume a signal without any "new physics" contamination, 
such as neutrino decays [31-33]. 

For this study, we will use this effective error on R to formulate the requirements to a 
neutrino telescope to be useful for terrestrial^ neutrino oscillation experiments, i.e., we 



""^Note that since the cross sections increase with energy, the actual event rates are determined by a 
balance between flux and cross sections. 

^Further on, we refer to the "terrestrial experiments" as the experiments with a terrestrial, i.e.. Earth- 



will show the results for different errors and read off the required precision. Of course, 
the effective error will not only depend on systematics, but also on the number of events 
expected. This signal height depends on the distance and energy release of the event, 
which means that it is not predictable until measured. Typically, several tens of events 
are expected for astrophysical neutrino sources, whereas of the order of hundred events 
can be expected for electron neutrino appearance at superbeams for large sin^ 2^13 ~ 0.1 
(see, e.g., Ref. [2]). This means that the actual statistical errors from astrophysical sources 
and superbeams (oc yN) may not be different by orders of magnitude, and a combination 
of data could make sense - a hypothesis to be tested in this work. We will typically use 
5%, 10%, and 20% errors on R in order to discuss the requirements to an astrophysical 
measurement in the following sections. For example, a precision of ~ 20% was found in 
Ref. [28] for a flux close to the Waxman-Bahcall bound [34, 35] measured by IceCube. This 
flux of E^ dNi, /dEy = 10^'^ GeVcm^ s~^ was assumed for one year at IceCube, which 
means that higher precisions may be expected for longer running times. Larger fluxes and 
therefore smaller errors on R could also be obtained for IceCube volume upgrades, optically 
thick sources, galactic sources, or sources not emitting protons at energies the Waxman- 
Bahcall bound is normalized to (see, e.g., Refs. [21,36]). We illustrate in Appendix 1X1 how 
to relate such errors on R to specific event rates in a neutrino telescope, where we also 
discuss principle challenges for backgrounds, systematics, and the different sources. 

Astrophysical source identification 

We assume that the type of the astrophysical source can be identified. The results of this 
study are then to be interpreted together with the assumption of that source. In order to 
identify the source, first of all note that the ratios R for the discussed sources are strongly 
separated even for different oscillation parameters, which means that a small error on R 
would instantly be a strong hint for the source determination {R ~ 0.26, ~ 0.66, and 
~ 0.5, respectively, for neutron beams, muon damped sources, and pion beams, as well as 
sin^ 26'i3 = 0). In addition, the energy dependence can be used for the source determination 
because the flavor composition may change as function of energy in a characteristic, source- 
dependent way, such as from a pion beam to a muon damped source [26]. Furthermore, 
associated signals in gamma or cosmic rays can be used for source identification (see, e.g., 
Ref. [18]). Different flavor ratios also help to identify the source (see, e.g., Ref. [37] even 
for arbitrary flavor compositions), and the Glashow resonance mentioned above may help 
as well. If an independent source identification is missing, only generic methods, such as 
in Ref. [19], can be used to infer on the neutrino oscillation parameters. In addition, we 
assume that there is little contamination from different production mechanisms producing 
different flavor ratios at the source. For example, neutrons can be generated by collisions of 
high energy protons on ambient photons and protons, or in the photo-dissociation of heavy 
nuclei. In the first case, the Ue from neutron decays would be negligible compared to the 
neutrino flux from pion decays, which means that we effectively have a pion beam source. 
In the second case, photo-dissociation of heavy nuclei produces a pure Ug flux, i.e., we have 
a neutron beam source. In either case, the obtained flavor composition at the source is 
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relatively clean. 

In summary, there might be three interesting sources of information for the neutrino os- 
cillation parameter measurements, which are very complementary, but all suffer from low 
precisions on a timescale of the coming ten years: Neutrino beams, reactor experiments, and 
astrophysical neutrino flux measurements. In this study, we demonstrate that the combi- 
nation of the complementary knowledge could allow for early information on the unknown 
neutrino oscillation parameters, while the knowledge and statistics from each individual 
source will be moderate or poor. Note that, in principle, other complementary information 
could be obtained from Oi//3/3-decay (for the mass hierarchy, see e.g. Ref. [38] and references 
therein), a galactic supernova explosion (see, e.g., Refs. [39-42]), or cosmology (see also 
Ref. [43] for a different combination of information). We will not discuss these sources at 
this place. 

This study is organized as follows: First, we illustrate and motivate the complementarity 
among superbeams, reactor experiments, and astrophysical sources in Sec. |21 Then we 
describe the simulation methods and assumptions used in Sec. El The results sections 
(Sec. mto Sec. |HI) are ordered by the timescale of relevance and level of technicality: In Sec. El 
we discuss possible early information on (5cp using reactor experiments and astrophysical 
sources only. In Sec. we illustrate the impact on the mass hierarchy measurements at 
superbeams and reactor experiments. Then in Sec. IHl we investigate CP violation and 
CP precision measurements. Furthermore, we demonstrate in Sec. [71 how and where in 
parameter space the potential to resolve the octant degeneracy could be improved. Finally, 
we discuss the impact of the measurement of all flavors in Sec. JHl which turns out to be 
very useful for pion beams. Note that we do not consider precision measurements of 9i2 by 
astrophysical fluxes, which have been discussed elsewhere [20]. 

2 Complementarity terrestrial-astrophysical 

In this section, we illustrate the complementarity between astrophysical sources and terres- 
trial neutrino oscillation experiments, i.e., neutrino beams and reactor experiments. Un- 
less noted otherwise, we use the following set of simulated neutrino oscillation parameters 
throughout this study [44-47]: Am^^ = 2.5 ■ IQ-^eV^ sin^ 2^23 = 1, Am^i = 8.2- IQ-^eV^ 
sin2 2ei2 = 0.83, and sin^ 2^13 = 0.1 somewhat below the CHOOZ bound [48], 5cp = 0, and 
a normal hierarchy. 

In order to study R = 0^/(0e-|-0r) for astrophysical sources, we assume that either the pro- 
duction or detection of the neutrinos is incoherent, or the mass eigenstates loose coherence 
while traveling because of the long distances (see, e.g., Ref. [49] for a formal treatment). 
Therefore, any oscillations and CP-violating effects average out in order to obtain 
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Note that the probabihties in these formulas are the averaged ones from Eq. (^J. Using the 
standard parameterization of the mixing matrix U [50] and our standard values of the other 
oscillation parameters, one can then calculate R as function of the oscillation parameters 
for the different sources. We expand R for the different astrophysical sources to first order 
in 6113 : 

^Neutron beam _ 0.26 + 0.30 ^13 COS (5cP , (4) 

^Muon damped _ 0.66 - 0.52 ^13 COS 5cP , (5) 

pPionbeam _ 0.50 - 0.14 ^13 COS (^CP • (6) 

Higher order terms in ^13 are relatively small (but not negligible). For neutrino beams, 
however, the (5cp-dependent terms in P^e are suppressed by the mass hierarchy, which means 
that the 6'^3-term is the leading term for large 9i^. At the first oscillation maximum and in 
vacuum, we find (see, e.g., Ref. [51]): 

P^e~2^2^±0.09^i3sin(5cP, (7) 

where the plus is for antineutrinos and the minus for neutrinos. Because most of the first- 
generation superbeams are operated close to the first oscillation maximum and have a very 
narrow beam spectrum, this approximation should be useful for qualitative discussions. 
Most importantly, neutrino beams are dependent on the CP-odd sin(5cp, whereas astro- 
physical sources are dependent on the CP-even cos(5cp. In addition, the amplitude of the 
signal for neutrino beams will be determined by 6'^3, whereas ^13 acts as a perturbation 
for astrophysical sources. The astrophysical information does not depend on neutrinos or 
antineutrinos, whereas neutrino beams usually exploit exactly this dependence by imple- 
menting neutrino and antineutrino channels. Finally, pion and muon-damped sources have 
a different sign in front of the cos 5cp-term than neutron beam sources, and therefore the 
roles of (5cp = and vr are exchanged. 

We show in Fig. ^ the exact dependence of the observables R (astrophysical neutrinos) and 
P^e (neutrino beams) on ^cp for different sources. While astrophysical sources have the 
largest modulation of the amplitude at ^cp = and vr because of the cos (5cp-dependence, 
neutrino beams are strongly influenced at (5cp = ±7r/2 because of the sin 5cp-dependence at 
the oscillation maximum. In addition, as discussed above, neutrino beams show a different 
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Figure 1 : Sources used in this study for which the signal depends on Sqp ■ We show the quantities R and 
Pfj.e {Pp.e), respectively, as function of (5cp for different values of sin 2^13. The shaded bars illustrate the size 
of the 5%, 10%, and 20% errors for the chosen central values (horizontal lines). For the terrestrial neutrino 
beam, we assume vacuum oscillations (or short enough baselines) and a measurement at the atmospheric 
oscillation maximum. 

behavior for tlie neutrino and antineutrino operation modes. However, note that the com- 
parison between neutrinos and antineutrinos can do only very httle close to 6cp = and 

TT. 

In order to illustrate the measurement precision of Sqp, we show possible error bars for 5%, 
10%, and 20% measurement errors as the shaded bars. From the projection of the curves 
onto these bars, we can immediately read off the required precisions for R and the relevant 
parameter regions in Sqp (for large sin^2^i3) for the different sources: For neutron beam 
sources, we may learn something about 6cp already for a 20% error on i? - provided that 
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Figure 2: Illustration of the synergy among superbeams, reactor experiments, and astrophysical fluxes 
(at the example of a muon damped source) in sin^ 26'i3-(5cp-space. Shown are the curves for constant 
rates (superbeam, reactor experiment) and constant R (astrophysical flux) going through the best-fit point 
sin^ 2^13 = 0.1 and (5cp = vr. 



5cp is close to or 7r. If considerably higher precisions can be achieved, other values of 5cp 
can be extracted as well. For muon damped sources, the required precision in R is similar. 
For pion sources, however, only precisions of the order of 5% can help. 

A different way to look at the synergy among superbeams, reactor experiments, and as- 
trophysical fluxes is illustrated in Fig. |21 In this figure, the curves for constant total rates 
(superbeam, reactor experiment) and constant R (astrophysical flux) going through the 
best-fit point sin^ 26'i3 = 0.1 and 5cp = vr are shown. For the superbeams, this means that 
without spectral information all the points on the respective curves are degenerate, which 
is a good approximation for narrow-band beams. The combination of neutrinos and an- 
tineutrinos still leaves a remaining degeneracy between and vr, which can, for large error 
bars, not be resolved by the reactor experiment either, because the degenerate solutions 
are still too close to each other. However, the astrophysical flux has a completely different 
dependence on sin^26'i3 and (5cp, which means that it could disentangle the two values for 
(5cp (0 and vr). Note that all the values on the reactor and astrophysical flux curves are 
exactly degenerate even when using energy information. 

As far as different measurements are concerned, we expect an impact of the astrophysical 
sources on CP precision measurements (especially for (5cp close to and tt), and for the mass 
hierarchy measurements at the superbeams because the sgn(Am3]^)-degeneracy is located at 
a different value of (fake) (5cp than the original solution. We have also tested the impact on 
sin^ 26'i3 exclusion and sin^ 2^13 discovery potentials, we we have not found any significant 
impact. The reason is that these measurements are dominated by the ratio of signal to 
background in the appearance channels of the beams, ie., the absolute event rate determines 



the performance. This event rate is, for neutrinos and a normal mass hierarchy, smallest 
close to ^cp ~ 7'"/2 (c/.. Fig. HJ lower right panel), i.e., ^cp ~ ^/2 (or 37r/2 for the inverted 
hierarchy) limits the performance. Since R for the astrophysical sources does not depend 
strongly on sin^ 26*13 close to 5q-p = tt or 37r/2, they are of little help (c/.. Fig. [H first three 
panels). Note that information on sin^ 2^13 cannot be expected from an astrophysical source 
alone either, because a different combination of 5cv and 6*23 can easily fake sin^ 26*13 = 0. 

Finally, note that R depends on ^23 in a different way than P^g for the beams (c/., Refs. [19, 
51]): While R is in some way proportional to ^23 by different combinations of powers of sin ^23 
and cos ^23; the main sensitivity for the neutrino beams comes from the term oc sin^ 6*23 ^13 
term in P^e- We will discuss this different dependence in Sec. [7| 

3 Simulation method 

In order to simulate the impact of R on the future terrestrial neutrino oscillation program, 
we use a modified version of the GLoBES software [52]. Let us define the input parameter 
vector used for GLoBES as 

'x = (6*12, 6*13, 6*23, 5cp, Am^i, Aml^, pi, . . . , pn) , 

where we simulate the combination of n terrestrial experiments with the average matter 
densities pi, i.e., we deal with 6 + n parameters. Note that the matter densities are treated 
as oscillation parameters in GLoBES in order to include correlations with matter density 
uncertainties. We implement the x^ as function of the simulated and fit parameter vectors 
as 

2 2 /— >■ — > \ I [R[^ Sim) " R{ X fit)) , , 

X =XTerr(a;sim, a;fit) + ^ , (8) 

where a^ = cr^' ■ R{ x sim) is the absolute error on R determined from an astrophysical 
source, and cr^' is the relative error which we will be using. In this equation, Xxerr i^ the x^ 
in GLoBES describing the terrestrial experiments and external input before marginalization 
over the oscillation parameters, and R is the flux ratio using the exact analytical expressions 
from Eq. (^. We marginalize x^ in Eq. (jHJ over any unused oscillation parameters as defined 
by the specific performance indicator (such as over all parameters except from 6cp for a 6cp 
precision measurement). Note that R does not depend on the mass squared differences and 
matter densities, which means that only four of the Q + n parameters determine the second 
term of Eq. (jH)). 

As far as the time scale for this study is concerned, we require substantial information 
from both a neutrino telescope and superbeam experiments. Therefore, we discuss neu- 
trino oscillation physics at the end of the first-generation superbeam era, which could be 
around ten years from now (until about 2015-2017). At that time, we assume that data 
from MINOS, Double Chooz, T2K, and NOj^A will be available for the part of the terres- 
trial experiments. In addition, there might be early large reactor experiments, such as the 
Double Chooz upgrade "Triple Chooz" [53], which we do not consider since they are less 
established yet. We use the MINOS simulation from Ref. [2] with a total luminosity of 
5yr X 3.710^°pot/yr and a 5.4 kt magnetized iron calorimeter [3] (the unit "pot/yr" refers 
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to "protons on target per year"). For Double Chooz, we use the simulation from Ref. [13] 
applied to Double Chooz, i.e., L = 1.05km and an integrated luminosity of 60 000 unoscil- 
lated events [11], corresponding to about three years of data taking with both near and far 
detectors. For T2K, we use the simulation from Ref. [54] based upon Ref. [4] with a total 
luminosity of 5yr x 1.0 10^^ pot /yr and a 22.5 kt water Cherenkov detector with neutrino 
running only. For NOz/A, we use the simulation from Ref. [55] updated to the numbers from 
Ref. [5] , i.e., L = 810 km operated 12 km off-axis with 3 yr x 6.5 10^° pot/yr for neutrinos and 

3 yr X 6.5 10^° pot/yr for antineutrinos and a 30 kt totally active scintillating detector (unless 
stated otherwise). In most cases, we will show the combination of these experiments, unless 
we use the reactor experiment alone for a somewhat shorter time scale. Note, however, that 
MINOS only has a marginal impact on the final results {i.e., the finally achieved MINOS 
luminosity hardly matters), and the main results will be mainly determined by accumulated 
statistics from the superbeams and Double Chooz, the mass hierarchy sensitivity of NOz/A, 
and the ability of Double Chooz to cleanly extract sin^2^i3. From the conceptual point of 
view, we therefore expect that the qualitative results will be hardly affected by the finally 
achieved luminosities as long as the latter experiments are present with sufficient statistics. 

Except from the information coming from the simulated experiments, we assume an external 
precision of 5% for each Atti^i and 612. This should approximately correspond to the 
precision at the time of the data analysis (see, e.g., Ref. [45]). In addition, we include a 5% 
uncertainty on the value of the baseline-averaged matter density, where the uncertainty takes 
into account matter density uncertainties as well as matter density profile effects [56,57]. 
Since a reactor experiment needs some information on the leading atmospheric parameters, 
we impose a 10% external error on Ar?7,|]^ for the analysis of a reactor experiment alone. 

4 Early knowledge on Sep from reactor experiments? 

The next generation of experiments providing some information on sin^ 2^13 should be re- 
actor experiments. In particular. Double Chooz could measure sin^ 2^13 already by about 
2011 or so [11] before the beams may have started running. Let us assume that sin^ 2^13 
is large, i.e., close to the current upper bound. In this case. Double Chooz could measure 
sin^ 2^13 relatively precisely (c/., upper right panel of Fig. |2I). However, the results will not 
depend on Sqp, which means that it is impossible to obtain information on 6cp from Double 
Chooz alone. In addition, it is difficult to obtain information on 6cp from an astrophysi- 
cal source alone as well after the other oscillation parameters have been marginalized over 
(c/., upper left panel of Fig. |2)). If, however, an astrophysical source is able to provide this 
information on a similar timescale as the reactor experiment, one will actually be able to 
learn something on 6cp already before the superbeams provide results. We show in Fig. El 
lower row, two examples for the combination of an astrophysical and Double Chooz signal. 
In both of these cases with a 10% precision measurement, a range of values of ^cp (arrows) 
could be excluded at the 90% confidence level. Of course, a CP violation measurement will 
not be possible because of the low precision and the low sensitivity of the astrophysical 
sources close to maximal CP violation. 

As far as the dependence on the true 6cp is concerned, Sep = and vr for neutron beam 
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Figure 3: Fit regions as function of sin^ 29i3 and 6cp for different experiments as given in the plot 
captions (the lower row is always in combination with Double Chooz). The simulated values are chosen as 
marked by the dots. The contours are shown for the la (black curves, dark regions) and 90% (gray curves, 
light regions) confidence level (1 d.o.f.). Dashed curves represent the results when the other (not shown) 
oscillation parameters are fixed, i.e., not marginalized over. The arrows in the lower row mark the ranges 
in Sep which can be excluded at the 90% confidence level. 

and muon damped sources yield similar qualitative results, For ^cp close to maximal CP 
violation, only muon damped and pion sources can provide some hint for excluded values 
of (5cp at the la confidence level for high precisions of R (about 5%). If the precision of 
the astrophysical measurement is only 20%, there will be some hints for some points at 
the la confidence level for neutron beam and muon damped sources. A 3a exclusion is 
only possible for a muon damped source if Sqp = 0. In this case, one can actually exclude 
6cp = TT at the 3a confidence level. 
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Figure 4: Fake solution Sep for the sgn( Am|i )-degeneracy as function of the true sin^ 2^13 for different 
values of the true Sep as given in the plot. Diamonds mark the points beyond which (to the right) the 
degeneracy can be resolved at the 90% confidence level. Figure for NOj/A (3 yr neutrinos+3 yr antineutrinos) 
and a normal hierarchy. 

5 Mass hierarchy determination 

The mass hierarchy determination using astrophysical flavor ratios alone will not be possible 
because the observables do not depend on the mass hierarchy for averaged oscillations. The 
first experiment from the next generation of beam experiments with sufficiently large matter 
effects to determine the mass hierarchy will be NOz/A, i.e., the mass hierarchy sensitivity 
will be dominated by NOz/A. However, the sgn(Am3]^)-degenerate solution of NOz/A, in 
general, appears at a different value of 6cp than the original solution, which means that 
the correlations with Sep and sin^ 2^13 highly affects the mass hierarchy sensitivity. We 
illustrate in Fig. E] the running of the fake solution in (5cp-space for this experiment for 
different simulated values of 6cp as function of the true sin^ 2^13. Note that the degeneracy 
can only be resolved on the right-hand sides of the diamonds of the curves (90% CL). One 
can identify two sets of fixed points in this graph: For small values of sin^2^i3, the fake 
solution is close to the original one (or vr — (5cp). For very large sin^ 2^13, the fake solution 
is close to an attractor point which is, for this beam, ^cp = 371/2 for the normal hierarchy 
and ^cp = ^/2 for the inverted hierarchy.^ As it is obvious from Fig. |3J a constraint on ^cp 

•^In order to understand this dependence and attractor points, the concept of bi-probability/bi-rate 
graphs is useful (c/., Refs. [8,58-60]). For small sin^ 2^13, the original and fake solution "pencils" overlap 
completely, and the fake solution Sqp is close to the original Scp or tt — ^cp (c/.. Fig. 2 in Ref. [60]). For very 
large sin 2^13, where the pencils are separated, the fake solution is close to an attractor point determined 
by the closest point of the fake solution ellipse to the original solution ellipse (for this beam, Sep = 37r/2 
for the normal hierarchy and Sqp = 7r/2 for the inverted hierarchy). Note that for T2K matter effects are 
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Figure 5: Sensitivity to the normal mass hierarchy as function of true sin^ 26*13 and true Sep (stacked 
to the "Fraction of <5cp") for MINOS, Double Chooz, T2K, and NO^A combined with an astrophysical 
neutrino source (as given in the plot captions) at the 2cr confidence level. The curves are for the following 
errors on the astrophysical flux ratio: no astrophysical flux observed (thick gray curves), 20% error on R 
(dashed curves), 10% error on R (thin black curves), and 5% error on R (thick black curves). 



would increase tlie confidence level of the degenerate solution if it was strongly running, i.e., 
very different from the original solution. Therefore, we expect an improved mass hierarchy 
sensitivity in combination with an astrophysical neutrino flux. 

In order to test this hypothesis, we show in Fig. the sensitivity to the normal mass 
hierarchy as function of true sin^ 26'i3 and true (5cp (stacked to the "Fraction of 5cp") for 
MINOS, Double Chooz, T2K, and NOz/A combined with an astrophysical neutrino source. 
The interpretation of this figure is as follows: For large sin^ 2^13, the terrestrial experiments 
alone will only be able to determine the mass hierarchy for about 50% of all possible values 
of 6cp which could be realized by nature. However, using, for instance, a neutrino beam 
source flux measured with a precision of 20% increases this fraction to 80% of all values 
of 6cp- Therefore, depending on chosen confidence level and precision of the astrophysical 
flux, the chance to discover the mass hierarchy will be improved from about half of all 
possible cases of 6cp to almost certain. Similar results can be expected from a muon 
damped source, while for a pion beam source the fraction of 6cp will be increased by up 
to 10%. Note that while MINOS, T2K, and Double Chooz combined do not have any 
relevant sensitivity to the mass hierarchy (c/., Ref. [2]), NOz/A alone could not achieve this 
impressive result in combination with the astrophysical sources either. Only the combination 
of these experiments will provide the information to disentangle sin^ 2^13, ^cp, and the mass 
hierarchy, as well as it will have sufficient statistics. In addition, we expect qualitatively 
similar results for the inverted hierarchy, where the role of Sep = 7r/2 and 3it/2 is exchanged 
for the superbeams, while the astrophysical sources still provide the relevant information 
close to (5cp = and n. 



too small to separate the two pencils sufficiently, i.e. 
shift the fake solution (still appearing at a low Ax^). 



the addition of the astrophysical source would only 
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Figure 6: Sensitivity to CP violation as function of true sin'^ 26'i3 and true 6cp (stacked to the "Fraction 
of 6cp") for MINOS, Double Chooz, T2K, and NOt^A combined with an astrophysical neutrino source (as 
given in the plot captions) at the 2(T confidence level (normal mass hierarchy assumed). The curves are for 
the following errors on the astrophysical flux ratio: no astrophysical flux observed (thick gray curves), 20% 
error on R (dashed curves), 10% error on R (thin black curves), and 5% error on R (thick black curves). 

6 Measuring Sep 

Learning about 6cp has two aspects: First, we want to measure leptonic CP violation. 
Second, we want to determine 6cp precisely. Note that these two options are not necessarily 
correlated: If, for example, 6cp is very close to (but not equal to) or n, we will certainly 
not be able to measure this small CP violation with the next generation of experiments. 
However, we may still be able to learn something about ^cp, such as we might exclude 
certain ranges. While superbeams are targeted towards CP violation, the sensitivity of 
the astrophysical sources should be best close to CP conservation. Therefore, we expect 
small effects for CP violation measurements, and larger effects for CP precision/exclusion 
measurements. 

We show in Fig.Elthe sensitivity to CP violation as function of true sin^ 2^13 and fraction of 
(5cp for MINOS, Double Chooz, T2K, and NOz/A combined with an astrophysical neutrino 
source. Note that the overall fraction of possible values for which this measurement could 
be done is rather small, i.e., between about 20% and 40% for large sin^ 26*13. The jump from 
20% to 40% comes from enough statistics to reduce the size of the sgn(Am3]^)-degenerate 
solution such that it does not overlap the CP conserving values anymore. Therefore, it is 
present depending on statistics and chosen confidence level (for instance, for NOz/A alone, 
it is not present). An astrophysical flux measurement could not increase the fraction of 6cp, 
for which CP violation is measurable, substantially (because it is most sensitive close to 
the CP conserving values), but it could increase the reach in sin^ 2^13 for about 20% of all 
values of 6cp. Note, however, that a similar effect could be achieved with slightly increased 
statistics. 

CP precision measurements are somewhat more complicated to describe, especially since 
the discussed experiments rather exclude some ranges of 6cp than they allow for a precision 
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Figure 7: The Ax^ as function of the fit value of Sqp for different simulated values of 5cp as given in 
the plot captions. The different curves correspond to no astrophysical information (thick gray curves), a 
20% measurement of R from a muon damped flux (thin curves), and a 5% measurement of R from a muon 
damped flux (thick black curves), all in combination with the terrestrial experiments. The dashed curves 
would respresent the solutions if the mass hierarchy degeneracy was resolved. The CP coverage is obtained 
as summed fit ranges at a specific confidence level, such as illustrated for the black curve and Ax^ = 9 in 
the left plot (gray bars). 

determination. For illustrative purposes, we show in Fig. [7| the Ax^ as function of the 
fit value of ^cp for two different simulated values. The different curves correspond to no 
astrophysical information (thick gray curves), a 20% measurement of R from a muon damped 
flux (thin curves), and a 5% measurement of R from a muon damped flux (thick black 
curves), all in combination with the terrestrial experiments. For (5cp = (left plot), the 
impact of the astrophysical source is obvious: the Ax^ for the solution at 5cp = vr is 
continuously lifted with increasing precision of R (arrow). Note that for a 5% precision 
measurement of i? a precision of 5cp is reached which allows to exclude maximal CP violation 
((5cp = 7r/2 and 37r/2) almost at the 3o" confidence level. For 5cv = 7r/2 (right plot), however, 
there is hardly any impact until high precisions of R are reached, when the Ax^ at ^cp = ti" 
is significantly increased. One can easily read off this effect from Fig. ^ (upper right): At 
(5cp = '7r/2 the error bars have to be about a factor of two smaller than at (5cp = or tt 
to exclude some ranges in ^cp- As for the CP violation measurement above, the effect at 
the 2a confidence level is small and the sgn(Am3]^)-degeneracy is already lifted above this 
confidence level at 5cp = vr (the difference between the dashed and solid curves describes 
the impact of this degeneracy). 

In order to study the CP exclusion/precision information as function of the simulated ^cp, 
we introduce the CP coverage [61]. The CP coverage is defined as range of fit values of ^cp 
which fit the chosen true value at a certain Ax^, and lies between (precise determination 
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Figure 8: CP patterns for MINOS, Double Chooz, T2K, and NOi^A combined for several external 
precisions on R from an astrophysical muon damped source. The CP patterns quantify the measurement 
of Sep (CP coverage) as a function of the true value of 5cp (provided by nature) for sin^ 26*13 = 0.1. The 
CP coverage is defined as range of fit values of 5cp which fit the chosen true value, and can be between 
(precise determination of Sep) and 360° (no information on Sep)- The thick curves correspond (from dark 
to light) to Ax^ = 9, 4, and 1 respectively, and the thin curves represent the results without taking into 
account the sgn( Am3i )-degeneracy. 

of 6cp) and 360° (no information on 5cp). It is illustrated in Fig. [7|as the gray bar in the 
left plot for the 5% precision measurement of R and Ax^ = 9. We use in Fig. |Hlthe "CP 
patterns" [60], which quantify the measurement of 6cp (CP coverage) as a function of the 
true value of 6cp provided by nature. This means that we compute Fig. [7|for all simulated 
values of 6cp, read off the CP coverage from each figure, and show this performance indicator 
as function of the simulated ^cp- Fig.|Hlshows these CP patterns for MINOS, Double Chooz, 
T2K, and NOi^A alone, as well as combined with a muon damped source for sin^ 2^13 = 0.1 
close to the CHOOZ bound [1]. Obviously, how precisely one can measure ^cf or what range 
of values one could exclude strongly depends on the chosen confidence level.^ However, it 
is obvious from this figure that astrophysical sources could especially help in the first two 
quadrants where the degeneracy problem is present (for a normal hierarchy, for the inverted 
hierarchy it is the third and fourth quadrants), and close to and tt. The CP patterns 
for a neutron beam source are qualitatively similar, whereas the pion beam source requires 
higher precisions. 

The structure of the patterns in Fig. |Hlis highly non-trivial. In particular, for the terrestrial 
experiments only, the precision is best close to ^cp = and tt at the la confidence level, 
and worst at these points at 3a. That is because the precision of Sqp depends very much 
on the size of the error bars. Going back to Fig. ^ (lower right) and studying the impact 
of the error bars (shaded bars, terrestrial neutrino beams), one can read off that for small 
errors, the precision will be best close to and vr as P^e leaves the simulated value very 
quickly. For large errors of the order of the modulation amplitude, however, the precision 



''Note that strictly speaking, one should refer to the Ax^ instead of the confidence level, since the 
connection between the shown Ax^'s and the Icr, 2a, and 3a CL for 1 d.o.f. holds only for quantities with 
Gaussian errors. 
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Figure 9: Impact of an astrophysical flux ratio measurement on the CP coverage for MINOS, Double 
Chooz, T2K, and NOi'A combined for the true Sqp = (left) and tt (right), as well as sin^ 2^13 = 0.1. The 
bars represent the la, 2a, and 3a measurements (Ax^ = 1, 4, 9) for different astrophysical sources as given 
in the plots. The right edges of the bars correspond to no astrophysical information, whereas the left edges 
correspond to a 5% measurement if the respective flux ratio. The vertical lines in the bars represent (from 
right to left) no astrophysical information, a 30% precision, a 20% precision, a 10% precision, and a 5% 
precision, where only some of the lines are labeled. 

will be worst close to and vr because at the values around the center easily all values 
of ^cp (for neutrinos and antineutrinos) can be covered (for more details, see Ref. [60]). 
Note that the size of the error bars depends on statistics as well as the chosen confidence 
level. Therefore, one expects a behavior strongly dependent on the confidence level for 
CP precision measurements (see, e.g., Ref. [61]). This argument can be translated to the 
astrophysical sources, where it will occur for and vr exchanged with ±7r/2. Therefore, we 
can understand this very complementary information from these two data sources. 

In order to compare the confidence level dependence and dependence on the astrophysical 
source more quantitatively, we show in Fig. El the impact of different astrophysical flux ratio 
measurements on the CP coverage for the true values of 6cp = (left) and vr (right). In 
this figure, the different vertical hues in the bars correspond to 5%, 10%, 20%, and 30% 
measurements of R (from left to right), as well as to no astrophysical measurement (right 
edges of bars). Selected precisions are marked with arrows and numbers. The results depend 
very much on the confidence level, value of 6cp, and flux used. However, one can read off 
this figure that even a 30% measurement of a neutron beam or muon damped flux can have 
a substantial impact (2cr, 6cp = 0). In addition, a 10% precision of these fluxes would allow 
for an exclusion of certain ranges at 3cr. Therefore, an astrophysical flux observation could 
be very useful to help restrict 6cp. 
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Figure 10: Illustration of the observables {R for astropliysical sources, total event rates for the beam) 
for the octant degeneracy resolutfon as function of sin^ 023 for sin^ 2^13 = (left) and sin^ 2^13 — 0.1 (right). 
The bands reflect the unknown value of ^cp- Note the scaling of the vertical axes depending on the source 
considered. The gray-shaded areas mark the 3a excluded region [47]. For NOz^A, we assume five years of 
neutrino running for this figure. 

7 Resolving the octant degeneracy 

As discussed in Ref. [19], the flux ratio R at neutrino telescopes has a distinctive dependence 
on 6^23 which is sensitive to the (6*23, vr/2 — 6'23)-degeneracy. Note, however, that without ad- 
ditional knowledge on the other mixing parameters, this information can only be extracted 
from R in very speciflc cases, but in a rather model-independent way [19]. We demonstrate 
in this section that the complementarity among superbeams, reactor experiments, and as- 
trophysical sources allows for an exclusion of the octant degeneracy for any substantial 
deviation from maximal mixing. For illustration, we show in Fig. ^J the observables {R for 
astrophysical sources, total event rates for the beam) for the octant degeneracy resolution 
as function of sin^ 6'23 for sin^ 2^13 = (left) and sin^ 2^13 = 0.1 (right). The bands reflect 
the unknown value of (5cp. In order to discuss the potential to resolve the octant degeneracy, 
compare the left branch of each source (sin^ 623 < 0.5) with the right branch (sin^ 6'23 > 0.5). 
If we assume that a reactor experiment determines sin^ 2^13 fairly well and the impact of 
the unknown 6cp is one of the main uncertainties, this picture should be quite accurate. 
For sin^ 26'i3 = 0, the situation is quite simple because, in this limit, 6cp is meaningless. 
The superbeam disappearance channel measures sin^ 26*23 very precisely, but it does not 
have information on the octant. The rates in the appearance channel are too low to imply 
any information, and reactor experiment is not affected by 623 (and therefore not shown). 
However, the flux ration R of the astrophysical sources is very different for sin^ ^23 < and 
sin^ 6*23 > 0, which means that a resolution of the octant degeneracy should be possible. 
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Figure 11: Exclusion of the octant degeneracy as function of sin^ 20i3 and sin^ 023 at the la (thin black 
curves), 2a (thin gray curves), and 3a (thick black curves) confidence level. The shaded areas shows the 
results from the terrestrial experiments only (MINOS, Double Chooz, T2K, and NOi/A combined), whereas 
the dashed curves represent the combination with a 20% precision measurement of a muon damped flux. 
The gray-shaded areas mark the 3a excluded region [47] . Figure for S(jp — and a normal mass hierarchy. 



However, for large sin^ 26'i3, the situation is different: The uncertainties on Sqp and sin^ 2^13 
(not shown) should affect the astrophysical sources. The beam, on the other hand, now has 
a substantial appearance rate with some information on the octant. We therefore conclude 
that the measurement for small sin^ 26*13 should be dominated by the astrophysical source, 
whereas the measurement for large sin^ 2^13 could be dominated by the beam (plus reactor 
experiment). The latter hypothesis needs to be quantified, because it it unclear from Fig.lTUl 
how the correlations with the other oscillation parameters affect the degeneracy. 

In order to demonstrate the resolution of the octant degeneracy quantitatively, we show 
in Fig. ^2 the exclusion of the octant degeneracy as function of sin^ 2^13 and sin^ 623. The 
shaded areas show the results from the terrestrial experiments only (MINOS, Double Chooz, 
T2K, and NOz/A combined), whereas the dashed curves show the combination with a 20% 
precision measurement of a muon damped flux as a representative example. Note that, for 
instance, NOz/A alone could not resolve the octant degeneracy because of the correlation 
with sin^ 26*13 and ^cp- However, it is the combination between NOz/A and Double Chooz 
which allows for the degeneracy resolution for large sin^ 2^13 (shaded contours). As expected, 
the terrestrial experiments alone cannot resolve the degeneracy for small sin^ 2^13, while 
the combination with a muon damped source (or similarly a neutron beam source) can. 
Note that this effect for small sin^ 26*13 is qualitatively comparable to the combination of 
atmospheric and long-baseline data using a megaton-size Water-Cherenkov detector [62,63], 
while it may be available at a shorter time scale. As far as the dependence on the simulated 
6cp is concerned, we observe similar small qualitative effects as in Ref. 
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Figure 12: Exclusion of the octant degeneracy (above upper and below lower curves) as function of 
the precision of R and sin^ ^23 at the 2a confidence level. The left plot corresponds to small (simulated) 
sin^ 2013 = 0, whereas the right plot is shown for (simulated) sin^ 2^13 = 0.1. The different curves are 
computed for the combination of different fluxes (neutron beam: thick curves; muon damped source: thin 
curves; pion beam: dashed curves) with the terrestrial experiments. In the right plot, the result from the 
terrestrial experiments alone is shown as well. The gray-shaded areas mark the 3a excluded region [47]. 
Figure for ^cp == and a normal mass hierarchy. 

Since the dependence on sin^ 26*13 is quite straightforward for all astrophysical sources, we 
focus in Fig. 1121 on the two physics cases sin^ 26'i3 = (small) and sin^ 26'i3 = 0.1 (large). 
The figure shows the dependence of the octant resolution on the relative error on R for 
different astrophysical sources. In the case sin^2^i3 = 0, where the terrestrial experiments 
alone do not provide any substantial information, a 10% to 20% precision of R from a 
neutron beam or muon damped flux would be very useful for a 10% deviation from maximal 
mixing. For a pion beam source, however, the precision has to be very high. Note that even 
for poor measurements of a neutron beam flux some information for significant deviations 
from maximal mixing (but still within the currently allowed 3a range) could be obtained. 
For large sin^ 26'i3, a 10% to 20% precision of R from a neutron beam or muon damped flux 
could improve the reach in 623 by a factor of two compared to the terrestrial experiments 
alone. However, beyond that the effect saturates quickly and the terrestrial experiments 
dominate the potential to resolve the degeneracy completely. Therefore, an astrophysical 
flux will be especially useful for small values of sin^ 26^13. Note that, compared to Ref. [62], 
we have not considered mixed degeneracies (wrong octant and wrong mass hierarchy) in this 
study, but we expect that one could exclude this degeneracy as well if the octant and wrong 
hierarchy alone can be excluded (similar to Ref. [62]). In addition, we have only considered 
the exclusion of the octant degeneracy, which is much more difficult than the exclusion of 
maximal mixing (which can be done with the disappearance channel at beams). However, 
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Figure 13: The ratio S = (t>e/4'r (representing the electromagnetic versus hadronic shower events) as 
function of <5cp for the astrophysical sources. The shaded bars illustrate the size of the 5%, 10%, and 20% 
errors for the chosen central values (horizontal lines). 

since the exclusion of the octant degeneracy is sufficient to exclude maximal mixing (but not 
necessary), we do know that there will be some sensitivity to deviations from maximal mixing 
as well. We expect this to be much weaker than the one coming from the disappearance 
channels of the beams via the sin^ 26'23-dependence independent of sin^ 2^13. 

8 What if one can actually measure all flavors? 

So far, we have assumed that we cannot discriminate Ve and Vr shower events and we have 
used R = 4)fj,/{4>e + 0t) as an observable. However, in principle, electromagnetic (from z/g) 
and hadronic showers (from Ur) can be distinguished by their muon content. In addition, for 
large enough energies, the tau track can be separated from the shower and may be identified 
by following double-bang and lollipop events (see Refs. [28,64] and references therein). Since 
we do not know the flux normalization, we can describe the neutrino oscillation physics by 
two independent quantities completely. We choose R, as before, and S = (j)e/4>r which could 
be extracted from the ratio of electromagnetic to hadronic shower events.^ We expect that 
this quantity might be determined with a lower precision than R, and it only becomes an 
observable for higher energies. This could make it an interesting additional observable for 
astrophysical pion beam sources. 

In order to discuss the qualitative dependence of the observable S on 6cp, we show it 
in Fig. ^1 for the different discussed sources. Note that we again illustrate the size of 
the relative errors as the shaded bars. Comparing Fig. ^J to Fig. ^ illustrates that for 
muon damped sources and pion beam sources the functions are decreasing in 6cp instead of 
increasing. Thus, a combination between 5* and R should compensate the difference from 
statistics between 6cf = and vr somewhat. Comparing the error bars in these figures 
Fig. m and Fig. ^1 there are, however, no major qualitative differences between S and R, 



^Note that neutral currents produce hadronic showers, too. In this case, we refer to the hadronic showers 
from Vt events only. 
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Figure 14: Relative improvement of the CP coverage for Sep = as function of the relative errors on 
S and R for the combination of the terrestrial experiments MINOS, Double Chooz, T2K, and NOz/A with 
an astrophysical pion beam flux {2a confidence level). 

which means that we do not expect additional synergistic information. In particular, for 
neutron beams and muon damped sources, the modulation of 5* is comparable to the one 
of R, which means that one would not expect major improvements except from increased 
statistics - especially given the fact that the expected precision of S may be weaker than 
the one of R due to lower efficiencies. The most interesting observation is probably the 
dependence of S on 6cp for the pion beam source, which is stronger than for R. Therefore, 
future precision measurements with high statistics (such as at IceCube volume upgrades) 
may find S a useful observable for standard pion sources. The same applies, in principle, 
to the dependence on 623 (not shown). We therefore focus on pion beams in this section. 

As a first example, let us study the dependence of the CP coverage on the relative errors of 
S and R from a pion beam source. Therefore, we show in Fig. E] the relative improvement 
of the CP coverage from the terrestrial experiments for ^cp = at the 2a confidence level. 
In this figure, the upper right corner (white area) corresponds to no improvement of the ter- 
restrial case and large errors on the astrophysical observables, whereas the lower left corner 
corresponds to very precise measurements of 5* and R. If we a assume a measurement of 
either S or R, we can read off this figure the required precision for a specific improvement. 
For example, for a 10% improvement of the CP coverage, a 27% precision of S (intersec- 
tion of 10% contour with right vertical axis) or a 17% precision of R (intersection of 10% 
contour with upper horizontal axis) are equally helpful. Obviously, the precision of S has 
to be substantially lower (by about a factor of two) than the precision of R to achieve the 
same result. For a 40% improvement, there is even a factor of three between the required 
precisions. One can also read off this figure that there are no real synergies between S and 
R: If both quantities are measured, one cannot expect an improvement beyond what is 
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Figure 15: Exclusion of the octant degeneracy as function of either the relative error on R or S, 
and sin'^ 6*23 at the 2a confidence level for sin^ 2^13 — 0. The curves represent the terrestrial experiments 
(MINOS, Double Chooz, T2K, and NOr^A combined) in combination with an astrophysical pion beam flux. 
The gray-shaded areas mark the 3a excluded region [47] . 

expected from the added statistics.^ 

As another example, we have indicated above that S could help for the resolution of the 
octant degeneracy. We illustrate its impact in Fig. ^1 where we show the exclusion of the 
octant degeneracy (same sign of A?7i|^) as function of the relative error on S* or i? at the 
2(7 confidence level. In this figure, sensitivity is given above the upper and below the lower 
curves. One can easily read off this figure that the precision of S has to be much lower 
than the one of R in order to do this measurement. Even a 30% error could eliminate some 
this degeneracy in a part of the currently allowed 3a allowed region of the parameter space, 
and a 10% error could eliminate the degeneracy for not too large deviations from maximal 
mixing. 

At the end it will be a matter of statistics which of the two observables S and R can be 
measured better. In principle, if each shower in the detector (used for R) was identified 
unambiguously, the precisions of S and R should be approximately comparable. That 
is because the neutrinos arrive approximately in the flavor ratio (1:1:1) for this type of 
source [27], which means that comparable relative errors of R and S should require similar 
event rates (somewhat depending on the efficiency ratio between track and shower events; 
see discussion in Appendix^. In practice, the precisions of 5* and R may be very different 
because of the required higher energies for the shower identification and backgrounds from 
mis-identification. Note, however, that the factor of two to three lower required precision 



^For Gaussian errors, the addition of two external measurements with an equal A^^ leads to about a 
factor of l/v2 smaller measurement error. 
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of S from physics corresponds to about a factor of four to ten lower event rate. Therefore, 
depending on the actual neutrino flux from the source (and the energies), S could be an 
interesting observable as well. 

9 Summary and discussion 

In summary, we have studied the requirements for an astrophysical flux measurement from 
a neutron beam source, muon damped source, and pion beam source at a neutrino tele- 
scope to be relevant for the neutrino oscillation program using terrestrial sources (reactors 
or accelerators). We have used the ratio R = 0^/(06 + <Pt) as our main observable, which 
can be extracted from the ratio between muon track and shower events in a neutrino tele- 
scope such as IceCube. This ratio might be measured at relatively low energies at IceCube, 
since the threshold for muon tracks is about 100 GeV, and for showers about 1 TeV.^ We 
have illustrated the complementarity of the observables R (astrophysical fluxes), Pgg (re- 
actor experiments), and P^e or P^g (superbeams) with respect to their 5cp-dependence: 
The observable of reactor experiments has no (5cp-dependence, the one of first-generation 
superbeams depends on sin^cp (at the oscillation maximum), and the one of astrophysical 
neutrino sources on cos 6cp- This complementarity has lead to a number of interesting obser- 
vations for large sin^2^i3: First, we have demonstrated that a 10% precision measurement 
from a neutron beam or muon damped source may lead to a first observation of 6cp together 
with Double Chooz, especially if the true 6cp is close to Sqp = or tt. This means that, 
depending on the time scale of the neutrino telescopes, Sqp may be restricted before the 
superbeams T2K and NOz/A start operation. Second, we have demonstrated that a 10% to 
20% precision of the ratio R = 0^/(0e + 4>t) coming from a neutron beam or muon damped 
source, or a 5% precision from a pion beam would be beneficial for mass hierarchy and 6cp 
measurements because of two factors: The cos (5cp-dependence of R improves the precision 
of Scp close to and tt, and the additional restriction on the parameter space prohibits the 
sgn(Am3]^)-degeneracy (mainly determined by NOz/A) to move in (5cp-space. For example, a 
moderately precise measurement of R from a neutron or muon damped source could enhance 
the chance to measure the mass hierarchy from about 50:50 to almost for sure. In addition, 
we have found that for small values of sin^2^i3, where the beams and reactor experiments 
alone do not have significant sensitivity, an astrophysical neutrino flux may exclude the 
octant degeneracy (for large sin^2^i3, reactor experiments and superbeams combined can 
already resolve this degeneracy). This application may provide constraints earlier than from 
atmospheric neutrinos and beams using a megaton Water Cherenkov detector [62]. 

Beyond the measurement of the muon track to shower ratio, we have discussed the impact of 
all-flavor identification. Therefore, we have introduced S = 4>e/4>T as an additional observ- 
able at higher energies (possibly using detector upgrades) which could be extracted from the 
electromagnetic to hadronic shower ratio. We have demonstrated that this observable could 
be especially interesting for pion beam sources, which might be the most likely sources to be 
observed at a neutrino telescope. The required relative precision of S has to be much weaker 



^However, it is difficult to differentiate electromagnetic (from i^e) from hadronic (from ly-r) showers at 
least below several PeV. Therefore, we have added the electromagnetic and hadronic shower events in R. 
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than on R for this type of source because of the stronger 5cp-dependence, where the role of 
S is similar to the one of R for (5cp, mass hierarchy, and octant degeneracy measurements. 

As far as the different sources are concerned, neutron beams and muon damped sources 
allow for the highest impact on the terrestrial program, which is, from the point of view of 
statistics, already present for several tens of events. However, neutrinos from neutron decays 
have a flux rapidly deceasing with energy and face the atmospheric neutrino background due 
to the steadiness of the source [23]. Muon damped fluxes, i.e., neutrinos from pion decays 
not followed by muon decays (which are absorbed before they can decay), occur in specific 
astrophysical models [25], and pion beams may turn into these sources at high energies [26]. 
Therefore, neutrinos from pion decays could be the astrophysical neutrino source with most 
statistics. Pion beams, however, have the highest requirement with respect to precision. 
There are several arguments in favor of these sources. First, as indicated above, shower 
(all-flavor) identification might help for pion beams. Second, because it may be the most 
likely sources, statistics from different sources and neutrino telescopes may accumulate. And 
third, once such fluxes are detected, fiducial volume upgrades of neutrino telescopes could 
easily increase statistics. Therefore, pion beams may not only have the highest requirements, 
but also the best prospects. 

With respect to time scales, IceCube [15] will be continuously deployed in the coming years, 
and other neutrino telescopes are in preparation as well [14, 16, 17]. Double Chooz may have 
accumulated significant data by 2011 or 2012, and the superbeams T2K and NOz/A may 
start data taking around that time. Beyond the coming ten years, other superbeam upgrades 
will follow, such as wide band beams [65] with a relatively wide energy spectrum. It is 
important to emphasize that astrophysical data can only be relevant for neutrino oscillation 
physics in the coming decade or at most fifteen years. For the discussed large values of 
sin^2^i3, mass hierarchy, and ^cp measurements will thereafter be certainly provided by 
superbeam upgrades having enough statistics (see, e.g., Ref. [66]), and the complementarity 
with respect to the cos 5cp-dependence will be provided by broad band superbeams, (3- 
Beams, superbeams operated off the oscillation maximum, or neutrino factories. For small 
values of sin^2^i3, atmospheric neutrinos and superbeam upgrades can exclude the octant 
degeneracy using megaton Water Cherenkov detectors [62,63]. 

We conclude that it is worth to study the ability of neutrino telescopes to extract the dis- 
cussed observables for neutrino oscillation physics, since this information may have a major 
impact on the neutrino oscillation program for the coming decade. Of course, the exact 
procedure and obtainable precision may require further research and the actual detection of 
a source. In addition, the interpretation of the data from all the discussed experiments will 
depend on the model of the astrophysical source. However, important hints for the planning 
of future experiments may be obtained early from the combination of a set of experiments 
with poor or moderate statistics each, but great synergistic potential. 
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A A note on statistics 

In this appendix, we give a simple example to illustrate the relationship between event 
rates and the error on R. We also discuss principle challenges for the different astrophysical 
neutrino sources. 

We assume that all simulated rates Tj and observed (fit) rates Oj are composed of signal 
events Si and background events Bi, i.e., 



th I f>th 
obs I riobs 



O, = Sr + Bf\ (9) 

where i corresponds to the muon tracks (? = 1) and (non-muon track) showers {i = 2). 
Now i?*^ = e~^ Sf^/S^^ and i?°^* = e^-*^ S'°^'^/S'|^* describe the simulated and observed ratios 
between muon tracks and showers. Note that these flux ratios have to be corrected by 
the ratio e of the efficiencies for muon track to shower detection e = ei/e2. Usually these 
efficiencies should be fairly well known. They depend on the fiducial volume (which can be 
different for different event types), cuts to reduce backgrounds, chosen energies, etc.. For 
the sake of simplicity, we will assume e = 1 in most cases, i.e., that muon tracks and showers 
are identified with the same efficiencies. For IceCube, one can easily derive e from Fig. 11 in 
Ref. [28] (see erratum) for the flux chosen therein. One finds e ~ 6 to 7 for a muon energy 
threshold of 100 GeV, and e ~ 3 for a muon energy threshold of 1 TeV (muon damped or 
pion flux; e ~ Muons/Showers ■ R~^). Therefore, we will use e = 7 for demonstration in 
some cases. Note that, in practice, the situation is more complicated and the mapping from 
incident to reconstructed neutrino energy needs to be described by the energy response 
function (usually a matrix, which describes the reconstruction at different energies). For 
example, some muon track events are partially contained (unless only the fully contained 
events are selected), which means that the muons are already loosing energy outside the 
detector. Thus, in general, e is a function of energy and therefore depends on the actual 
flux. In addition, since the muon range depends on energy (see, e.g.. Fig. 1 in Ref. [67]), the 
effective fiducial volume for partially contained events is energy dependent, as well as the 
event rates somewhat depend on the flavor composition. Here we assume that the energy 
response of the detector is known, and that e represents the integrated efficiency ratio for 
the flux observed.^ 



^In practice, the efficiencies for electromagnetic and hadronic shower detection are somewhat different, 
too, which means that the actual observable is R' = (t)^l{k4)e + </'r)- From Fig. 11 in Ref. [28] one can 
extract k ~ 0.7 for the flux and cuts chosen therein. For this value, the relative error on R' translates into 
the relative error on R by multiplication with a factor of ^ 1.1 to 1.3 for the discussed sources, whereas the 
qualitative dependence of R' on sin 2^13 and 5qp stays the same. Since this relative efficiency depends on 
the instrument and analysis, we do not use it in this study. 
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Let us introduce some simple systematics, too: 

Oi = (l + a) + (l + 6)5i, 

o^ = t:^ + (i + ^)^^' (10) 

where a, b and c are systematics auxiliary variables to be minimized over, and Bi = Bf^ 
are the predicted background events for muon tracks and showers. Thus, a describes the 
unknown normalization of Si (arbitrary, therefore we have not included it in the formula), 
which is, however, fully correlated between muon track and shower events. Furthermore, 
h and c describe the background normalization for muon tracks and showers, respectively. 
Now we can define a very simple Gaussian x^ as 

' {T,-0,f , fh\\ fc^' 



i=l 



,^ = 53iii^^ + ^ + i , (11) 



where ab and ac are the muon track and shower background normalization errors, respec- 
tively, and the auxiliary variable a remains unconstrained. Now we can compute the error 
on R°^'^ for certain assumptions for the event rates Ti and T2 (related by -R*^) and for 
systematics by marginalization over the systematics auxiliary variables a, b, and c. 

As far as backgrounds are concerned, these can be very different depending on the type 
of the source. Let us focus on the atmospheric neutrino background in this example. The 
atmospheric neutrino background makes several ten thousands of i/^ events per year for 
El, > 100 GeV, several thousands per year for Ey > 1 TeV, and about one hundred per year 
for Ey > 10 TeV (see, e.g., Ref. [68]). Therefore, it clearly is a problem at low energies. 
Now there are different principle strategies to reduce this background, such as: 

1. Good angular resolution 

2. Short timescale of event 

3. Cut low energy events 

For IceCube, the good angular resolution works for muon track events (~ 0.7°), but not 
for shower events (~ 25°) [28]. However, other neutrino telescopes may have a better 
angular resolution. The second approach helps for all short-timescale events. For instance, 
for gamma ray bursts, the atmospheric neutrino background should be negligible because 
of the short duration. The third approach can be useful for any source producing very 
high-energetic neutrinos. In summary, these strategies could be useful for many types of 
sources, such as muon damped fluxes. However, the observation of R for the neutron beams 
is particularly challenging because of the steadiness of the source and the flux dominating 
where the atmospheric neutrino background is large [23]. Therefore, we conclude that the 
specific systematics needs further study and clearly lies beyond the scope of this work. 

We show in Table ^ the event rates for muon tracks and showers for different hypotheses 
and systematics assumptions. Note that, as opposed to what is assumed for the rest of 
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this work, the errors are often not symmetric. As one can read off this table, backgrounds 
become important as soon as the background rates are of the order of the signal rates. For 
background free measurements and e = 1, we find, depending on the source, a minimum 
of about 30 to 50 events necessary for a 20% precision measurement, 100 to 150 events for 
a 10% precision measurement, and about 500 events for a 5% precision measurement. In 
this case, the statistics is determined by the muon track event rate, whereas for larger e the 
shower event rate becomes more important. 
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+0.20 
-0.17 


+0.47 
-0.26 


+0.33 
-0.23 




100 


152 


22 


+0.14 
-0.12 


+0.30 
-0.20 


+0.18 
-0.15 




250 


379 


54 


+0.08 
-0.08 


+0.17 
-0.13 


+0.09 
-0.09 


+0.40 
-0.23 


500 


758 


108 


+0.06 
-0.06 


+0.12 
-0.10 


+0.06 
-0.06 


+0.18 
-0.13 


Pion beam (R*'' = 0.5) 








10 


20 


3 


+0.45 
-0.35 








20 


40 


6 


+0.30 
-0.25 








30 


60 


9 


+0.24 
-0.21 


+0.56 
-0.30 


+0.44 
-0.29 




50 


100 


14 


+0.18 
-0.16 


+0.40 
-0.25 


+0.27 
-0.21 




100 


200 


29 


+0.13 
-0.12 


+0.26 
-0.18 


+0.16 
-0.14 




250 


500 


71 


+0.08 
-0.08 


+0.15 
-0.12 


+0.09 
-0.08 


+0.29 
-0.19 


500 


1000 


143 


+0.06 
-0.05 


+0.10 
-0.09 


+0.06 
-0.06 


+0.14 
-0.11 



Table 1: Some event rates for muon tracks, the corresponding predicted shower rates for different 
efficiency ratios e, and relative errors Ai?/i?*'* (Icr) for specific assumptions for backgrounds. The columns 
"No BG" and "e — 7" are calculated background-free for e = 1 and e — 7, respectively. The columns 
"10/100 BG" and "100/1000 BG" are calculated with backgrounds for e = 1, where the numbers refer to 
muon track background events/shower background events. In addition, a 10% background normalization 
error uncorrelated between muon tracks and showers is assumed for these columns (ab = cFc = 0.1). 
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